Renormalization Group study of the Sliding Luttinger liquids 
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We derive the RG-flow equations of the sliding Luttinger liquid perturbed by charge- density- 
wave (CDW) and superconducting (SC) operators. Using them we study the phase diagram of 
an array of XXZ spin chains coupled by Ising terms. In the weak coupling regime we find a new 
class of non-gaussian and unstable fixed points whose existence is due to a balance between the 
CDW, SC and sliding Luttinger couplings. 
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The Luttinger liquid (LL) is the key concept to de- 
scribe interacting electrons in one dimension. The collec- 
tive nature of the excitations, the spin-charge separation 
and the power law behaviour of correlators with anoma- 
lous exponents are some of its distintive feautures, in 
contrast to the Fermi liquid theory pj] . Much of the the- 
oretical activity in the last years has been devoted to the 
search of non-Fermi liquid theories in dimensions higher 
than one, specially in 2D, due to its possible connection 
with the high-T c superconductors and another strongly 
correlated systems. 

A natural path to realize a non- Fermi liquid has been 
to couple arrays of ID Luttinger liquids forming ladders 
and 2D planes. However the general consensus has been, 
until recently, that 2D arrays of LL's are unstable to the 
formation of crystal, superconducting or 2D Fermi liquid 
states . An alternative to these "no-go theorems" has 
been proposed lately using the concept of sliding Lut- 
tinger liquid (SLL), also called smectic non- Fermi liquid 

111- 

These works were partially motivated by the Ander- 
son's proposal of confinement of excitations in the Lut- 
tinger liquids || and have a classical analogue in the 
stacks of coupled 2D XY models coupled by gradient in- 
teractions @. The SLL model may also be relevant to 
the stripe phases of the Quantum Hall effect and the 
cuprates. 

The sliding Luttinger liquid is the fixed point of a gaus- 
sian Hamiltonian which treats on equal footing the in- 
dividual Luttinger Hamiltonians of the stripes and the 
density-density and current-current interstripe interac- 
tions ||, Q|. Using bosonization techniques one can re- 
gard the SLL as a set of decoupled LL's characterized by 
a sound velocity v(q±) and a Luttinger coupling K(q±_) 
which depend on the tranverse momentum q± across the 
stripes. The gaussian nature of the SLL fixed point allows 
a simple derivation of the scaling dimensions of the single 
particle (SP), charge density (CDW) and superconduct- 



ing (SC) order parameters, in terms of the SLL function 
K(q±) H f|. These scaling dimensions have been used 
to study the stability of the SLL under various perturba- 
tions, finding rich phase diagrams where the SLL phase 
survives in the vicinity of CDW, SC and Fermi liquid 
phases || ||, || . The perturbative RG analysis performed 
in the latter references take into account the running of 
the coupling constants of the perturbations to first order, 
while the Luttinger functions K(q±) and v(q±) stay con- 
stant. It is however well known that in certain situations, 
as in the presence of marginal perturbations, one has to 
consider in addition the renormalization of the Luttinger 
parameters in a RG a la Kosterlitz-Thouless (KT) 0. 
This is for example the case of the Hubbard model at half 
filling where the Umklapp operator becomes marginally 
relevant for a repulsive Hubbard constant, leading to a 
charge gap in the low energy spectrum. 

The aim of this Letter is to derive the RG equations 
of the Sliding Luttinger liquid for the coupling constants 
and the SLL functions to second order in the couplings, 
and study some of their consequences in a model con- 
sisting of arrays of spinless ID Hamiltonians coupled by 
density-density interactions. The latter model has been 
treated in the past with bosonization [^) , mean field E3] 
and variational methods |ll| , which predict the existence 
of large regions in parameter space where the phase is ei- 
ther CDW-like or XY-like, corresponding respectively to 
the smectic crystal and the smectic metallic phases of 
references [|[ Q. 

On more general grounds we also analyse the stabil- 
ity of the spinless SLL under marginal CDW and SC 
perturbations, reaching the conclusion that when one of 
these perturbations is irrelevant the other one is relevant 
and hence the SSL is an unstable fixed point. We also 
find new unstable fixed points if two conditions are sat- 
isfied: i) when the CDW coupling constants are minus 
the superconducting ones, which guarantees the freezing 
of the SLL parameters, and ii) when the running of the 
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CDW and SC parameters to first order, which is given 
by their scaling dimension, is cancelled out by the second 
order term. The mechanism involved here is similar to 
the one giving rise to the Wilson-Fisher fixed point in 
dimensions d < 4 Jl2] ]. From a more mathematical view- 
point the manifold of these non gaussian fixed points are 
closely related to the target manifolds that appear in 5 
dimensional supergravity theories |D| . 

Let us consider an array of N spinless Luttinger stripes 
with phase fields for the density fluctuations <j) a (a = 
1 , . . . , N) and euclidean Lagrangian 
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N Tf 



-{dt^af +V (d x <j) a ) 2 



(1) 



where Kq is the inverse of the standard Luttinger param- 
eter ( Kq > 1 for repulsion) and vq is the sound velocity 
which we scale to 1. The charge density fluctuations j'q 
and the charge currents jf are given by the bosonization 
equation = ^e fiu d u 4>a- Hence the density-density and 
current-current interactions among the stripes are also 
quadratic in the derivatives of the bosonic fields and, to- 
gether with (P, define the sliding Luttinger Lagrangian 



where K(q±) 

tinger parameter and v(q± 



y/K J (q±)KP(q±) is the (inverse) Lut- 
y/Kp(q ± )/KJ(q ± ) is the 



sound velocity of the q±~ mode. In the dual variables 
6 q± the Lagrangian has the same form as eq. (|J) with the 
replacement K(q±) — ► 1/K(q±). The scaling dimensions 
of the various operators, given in references j| @, can be 
easily derived by observing that a generic vertex opera- 
tor Vp = exp(i^2 a a (f> a ), factorizes into the product of 
vertex operators on the q±— modes as follows 
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E „ 0a (cos^a)^ 1 ' + s m(q ± a)<t>W ) 
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where (frq 1 ^ are real scalar fields. The scaling dimen- 
sion of Vg is given simply by the sum of the scaling 
dimensions of all the individual q±— modes, namely 
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where K 
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K are N x N matrices whose 



off-diagonal elements are given by the interstripe current- 
current and density-density interactions. The SLL model 
can alternatively be formulated in the dual variables 9 ai 
which arc the phase fields of the superconducting fluctu- 
ations. The SLL Lagrangian (Q) becomes || [|, || 



1 N 



„b=l 



(3) 



where Kj are the inverse matrices of K J ' P . Both 
cqs.(|^) and (|^) exhibit the smectic or sliding symmetries 
4>a — > 4>a + a a and 6 a — > 9 a + a , where a a and f3 a 
are constants, which prevents the lock in of the charge- 
density-wave and superconducting order parameters of 
the individual stripes Q]. Assuming periodic bound- 
ary conditions across the stripes and translational invari- 
ance along them, one can perform the Fourier transform: 



<t>a = -^E q± e^ a <f> q± , K J >?(q ± ) = J2 q± e^ a K^ +a 
in order to bring the SLL Lagrangian (||) into the form 
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\d t <f> q± \ 2 +v(q ± )\d x cb q± \ 2 



(4) 



where J — Similary, a vertex operator in the 

dual variables, i.e. Vi — exp(i J2 a f3 a 9 a ), has a scaling 
dimension Ajg given by the formula (||) with K(q±) re- 
placed by its inverse. The interaction Lagrangian is given 
by the pair hopping (SC) and particle-hole (CDW) oper- 
ators§,@, 



/ (2toq) 2 ^ [ 9 cdw cos/3(0 Q - <j) b ) (7) 



+g%' c cosf3(d a 



where ao is the lattice spacing and (3 — \/2tt for the 
charge modes of spin gapped systems j|, or (3 = y/iir 
for spinless fermions Q. In the latter references the sta- 
bility of the perturbed SLL was studied in terms of the 
relevance or irrelevance of the CDW and SC operators 
([?]), given by their scaling dimensions. There are how- 
ever cases where one has to consider the renormalization 
of the functions K(q±) and v(q±), for example when the 
CDW and SC operators become marginal. The latter 
situation actually arises in the study of the coupled XXZ 
spin chain Hamiltonians via Ising and spin-pair-flipping 
terms § 0, 
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(8) 

+ J'z^i,a^i,a+l + J'xY{^i^a^i+l,a^i : a+l^i+l,a+l + h.C.)] 
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This model can be Jordan- Wigner transformed into a 
spinless fermion Hamiltonian which, after bosonization, 
becomes at half filling B , 



N 

H = ^2 dx 



+ -7— (d x (j) a ) 
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2J' z a 
(2™ ) 2 



(cos\/47r(<?(> a + 4> a+1 ) - cos\/47r(0 a - a +i)) 



where u = Jao(l + 2J Z / (ttJ) and JTo = 1 + 2 J z /(tt J) for 
| Jj | < < J. The gaussian terms in ^ yield a SLL model 
(0) with 



where l\ a x b are the scaling dimensions of the CDW and 
SC operators ( eq.(||)) 



a, b 

^SC/CDW 



^-J (l-cos^a-fc))^ 1 ^) (13) 



and the functions and are defined as, 



t-iO \ ^ / a,b\2 T7ia,b \ ^ a,c c.b 

F x = l^V9x ) > F x =2-^9x 9x 



(14) 



As an application of these eqs. we shall consider the case 
when K(q±) = 1 + k + J2 n >o k n cos(q±n) is close to 1 
and (3 = -\/47r, which corresponds to the spinless fermion 
model (|[ [)]). The linearization of eqs.(|l2|) yield 



6 a , b Ko, K f B 
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(10) 



where the time variable and the exchange couplings are 
measured in units of u and J respectively. In the weak 
coupling regime \J Z \ << J, the SLL function K(q±) is 
close to 1. Consequently the intra-chain umklapp cou- 
plings cos-v/Tfhn^a have dimension d ~ 4 and hence can 
be neglected, as it is the case of decoupled single chains. 
On the other hand the CDW, SC and Umklapp inter- 
chain couplings are marginal and one has to consider 
their running together with that of the SLL functions. 
To simplify matters we shall neglect in what follows the 
Umklapp term, which is anyway absent away from half- 
filling. The couplings constants in (0) , which correspond 
to @, are given by 



a, b 

9cdw 



- 2 -^|a-6|,l! 



9SC = %J'xY<i\a-b\,l (11) 



Eqs.(||) and (|?|) define a multicomponent sine-Gordon 
model which can be renormalized using operator prod- 
uct expansion (OPE) techniques 12 . We shall further 



assume that the sound velocity of the modes is constant 
and set it to 1, so that the hard disk regularization will 
be common to all the modes. The RG-flow equations for 
K «,b = K J J and 9cDW/sc,n are given by ( up to second 
order in the g's) 



dK n 



ff 2 
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^ =- g^idCD W,n~ 9sc, n)> (n>0) (15) 

dgcD\v,n , ^ 1 dJV(gcDw) 
- (2k Q - k n )gcDW,n 



ds 
dgsc,n 



-(2fco - 



47r dgcDW, r 
1 dN{g sc ) 



ds ' "" 4^ dgsc,r 

where N{g) is the cubic polynomial 



A%) 



(16) 



n,m>0 



which encodes the OPE of the CDW and SC operators 
appearing in (^). Notice the invariance of (|l6|) under 
the change g n — » (— l) n g n . From ( |l5| ) it follows the RG 
conservation of K (q± = 0) = 1 + ^„> k n . Moreover the 
"Minkowskian" norm /z, defined as, 



n>0 

satisfies the RG eq. 

d/i 6 
= (4tt) 2 



(4tt) 5 



^ (9cDW,n + 9 2 sc,n) (17) 



n>0 



W(gcDw) +N(gsc)} 



(18) 



Eqs.(|l5|) define an RG-flow in an infinite dimensional 
space formed by all the SLL, CDW and SC parameters. 
To extract some information from them we have to do 
further approximations. Let us first consider the simplest 
case where only fco and gcDW,i are non vanishing. A 
truncation of eqs.(|l5|), written in the variables x = —2kg 
and y = gcD\vs/{^), gives 
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dx 2 dy 
ds ^ ' ds 



-xy 



(19) 



which are the well known Kosterlitz-Thouless RG eqs. 
of the XY model @, [l2|. In this case the norm fi ([l7| ) 
is proportional to x 2 — y 2 , and it is conserved since the 
cubic term vanishes. The RG trayectories are hyperbolas 
of constant /i which end up at infinity provided \y\ > 
x. The latter inequality translates in the XXZ+Ising 
model (||) into the condition | J' z \ > — AJ Z , which coincides 
with the one obtained in reference || for the existence of 
an AF-Ising phase. In the weak coupling regime of the 
XY model, which corresponds to the region \y\ < x, the 
parameter y flows to zero and x flows to the line x > of 
SLL stable fixed points. The inclusion of the interstripe 
Luttinger parameter ki modifies the boundaries of the 
strong coupling region. Indeed performing the change of 
variables x = — 2ko + fci and y = V6 gcD w,i/(47r) one 
obtains from ( |i"5j ) the same KT eqs.(|l9|), and the strong 
coupling regime is given now by J' z > —17.798 J z ( if J' z > 
0) and - J' z > -1.798J Z ( if J' z < 0). In all these cases a 
necessary condition to achieve a XY or SLL phase is to 
have a ferromagnetic intrachain coupling J z < 0, which 
agrees with the results of Q] concerning the proximity 
of the stable SLL to the isotropic ferromagnetic point 
where the boson stifness Kq vanishes. To consider the 
effect on the higher order couplings we make the change 
of variables x = — 2fco and y n = gcDW,n/ {^), obtaining 
the eqs. 



dx 
ds 
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Erf 



dy n 
ds 



-xy„ 



1 dN(y) 

2 dy n 



(20) 



If the initial conditions are set to x(0) — x, yi(0) 



This equation has non trivial solutions for g n as func- 
tions of k n . The fixed points we have analized turn out 
to be unstable. The non-gaussian nature of these fixed 
points recall the Wilson-Fisher fixed point, with 2fco — k n 
playing the role of e = 4 — D. From a completely different 
viewpoint it is worth to mention that the mathematical 
structure underlying these fixed points, is very close to 
the one that appears in 5D-Supergravity theories [ O , 
where the manifold of the scalar fields, that come in the 
Maxwell multiplets, is defined in terms of a cubic norm 



Af similar to the one introduced in (16). This connection 



may help to study the properties of the new solutions. 
An interesting problem is to investigate the existence of 
stable non-gaussian SLL fixed points when K(q±) does 
not lie near 1. This can be done starting from eqs.jl^) 
rather than from their linearized version (|l5|). Another 
issue is to include the effect of non constant v(q±). The 
generalization of our results to SLL's with charge and 
spin degrees of freedom is rather straightforward. Here 
too we expect the appearance of novel non-gaussian fixed 
points. 
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